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Two-dimensional graphene, carbon nanotubes and graphene nanoribbons represent a novel class 
of low dimensional materials that could serve as building blocks for future carbon-based nanoelec- 
tronics. Although these systems share a similar underlying electronic structure, whose exact details 
depend on confinement effects, crucial differences emerge when disorder comes into play. In this 
short review, we consider the transport properties of these materials, with particular emphasis to 
the case of graphene nanoribbons. After summarizing the electronic and transport properties of 
defect-free systems, we focus on the effects of a model disorder potential (Anderson- type), and il- 
lustrate how transport properties are sensitive to the underlying symmetry. We provide analytical 
expressions for the elastic mean free path of carbon nanotubes and graphene nanoribbons, and 
discuss the onset of weak and strong localization regimes, which are genuinely dependent on the 
transport dimensionality. We also consider the effects of edge disorder and roughness for graphene 
nanoribbons in relation to their armchair or zigzag orientation. 
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I. INTRODUCTION 

Since the discovery of fullerenes (Cen) jCarbon-based materials have been the subject of intense research, which led to 
the discovery of carbon nanotubes [l|, M, S 3 and the fabrication of individual one-atom thick graphene layers @ . This 
opens unprecedented avenues for the investigation of quantum transport in low dimensional ID and 2D systems, as 
well as attracting the interest of industries, given the potential for innovative applications. Carbon nanotubes science 
is now a mature field of research, whose theoretical foundations have been reviewed in many textbooks (see [3, 0] 
and references therein). More recently, the possibility to single out a graphene plane, either through an exfoliation 
process [H| , or by means of an epitaxial growth mechanism , has allowed the whole community of mesoscopic physics 
to revisit basic foundations of quantum transport in 2D systems, such as quantum Hall effects or weak localization 
phenomena 0, H, 0, [13, [HI, [13, ES, [l3| , in a material with remarkable electronic properties Additionally, the routes 
for an alternative carbon-based nanoelectronics are actively investigated ranging from device optimization to graphene 
integration at the wafer scale 0, [l^, . [T^. [T§. [isl. [20| . However, since truly two-dimensional graphene is a zero-gap 
semiconductor, its use as an active electronic device, such as a field effect transistor, requires the reduction of its lateral 
size to benefit from induced quantum confinement effects. Energy band-gap engineering has been demonstrated with 
both types of materials, carbon nanotubes (CNTs) and graphene nanoribbons (GNRs) [U, [IS]- In the case of 
GNRs-based field effect transistor, this results in an increase of performances with downscaling the ribbon width 
from several tens of nanometers to 2 nm [l^. [itI. [isl. [l9l. [20l. [23| . In contrast to carbon nanotubes, whose integration 
in operating devices at a large scale remains a tremendously complicated challenge, GNRs present the potential 
to be fabricated and massively integrated in complex architectures, thanks to the use of conventional lithographic 
techniques. 

In this work, we discuss the transport properties of carbon-based low dimensional materials by focusing on the 
effect of a case of disorder potential, namely the Anderson type model [131 ■ This model has been widely employed in 
studying the scaling theory of localization, and through analytical and numerical results, it allows us to overview the 
behavior of the basic transport length scales (elastic mean free path and localization lengths) and associated transport 
regimes in systems with different symmetries and transport dimensionality. In the case of graphene nanoribbons, the 
presence of the edges exposes the system to further sources of disorder. In order to account for this peculiarity, we 
also consider the effects of impurities and roughness at the edges of the GNRs on the transport properties. 

Section In] is devoted to an overview of the electronic structure and quantum transport properties of 2D-graphene, 
quasi-one-dimensional carbon nanotubes and graphene nanoribbons in the absence of disorder. We start from a simple 
nearest-neighbors (n-n) tight-binding model, and highlight differences and similarities between these different carbon 
structures due to their different dimensions and lack or presence of edges, as these are the underlying factors that 
determine the specific response of their transport properties to disorder. We include a discussion about how the 
electronic structure of narrow graphene nanoribbons is modified when performing ab initio calculations and how to 
mimic this behavior by means of a third nearest-neighbors tight-binding model. In the case of nanoribbons, we also 
elaborate on the dependence of spatial symmetry of the system on the armchair or zigzag edge orientation and the 
number of carbon atoms in the primitive cell, since in certain cases, this has a deep impact on the transport features. 
At the end of the section, two interesting examples of transport in clean nanoribbons are also analyzed. The first one 
is the so called "pseudodiffusive" transport regime, which presents curious and intriguing analogies with the diffusive 
transport regime of disordered metals. The second one explores the effects of high magnetic fields on the energy 
bands, the transport properties and on the spatial distribution of electrical currents fiowing through the ribbons. 

In section IIIIl we consider the effect of disorder on the density of states and on the transport properties of the 
carbon nanotubes and nanoribbons. We focus on the analysis of typical length scales such as the mean free path and 
the localization length, taking an Anderson-type disorder as a reference model. In spite of its simplicity, this type 
of disorder allows us to point out the differences between 2D graphene and other graphene-derived materials. More 
realistic types of disorder would lead to more complicated effects that might mask the essential dependence of the 
transport properties on the dimensionality and symmetries of the investigated systems. A subsection is also devoted 
to the edge disorder (roughness and Anderson-type) in graphene nanoribbons. The effect of this kind of disorder 
on the transport properties very critically depends on the armchair or zigzag edge orientation of the ribbons. In 
the literature, several different models of edge disorder are considered and investigated with not always unanimous 
conclusions. In the last part of the section, we discuss the charge mobility in graphene based systems subjected to 
different types of disorder, since such quantity is key to assess the performances of device performances. 

Section Uv] concludes the report by stressing the most relevant features. 
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II. ELECTRONIC PROPERTIES AND TRANSPORT IN DEFECT-FREE SYSTEMS 

In most cases, the energetics of graphene based materials can be conveniently described as a first approach by a 
TT-effective electron model Indeed, for the honeycomb geometry, a bonds are formed by three out of four valence 
electrons (sp^-hybridization), whereas the remaining single electron per carbon atom occupies a 2pz orbital whose 
hybridization with first neighbors generate tt and tt* bands. A restriction of the Hamiltonian to a single Pz orbital per 
carbon atom provides thus an effective energetic model that well describes the electronic properties close to the charge 
neutrality point (CNP). In the case of graphene nanoribbons, we assume that the a bonds at edges are saturated with 
H atoms, thus eliminating dangling bonds. 



A. 2D graphene band structure 



CNP locates the Fermi level for 2D undoped graphene and metallic nanotubes and graphene ribbons, and the 
mid-gap of semiconducting CNTs and GNRs. Defining 70 as the integral overlap between Pz orbitals, the electronic 
properties of 2D graphene are straightforwardly derived by diagonalizing a 2 x 2 matrix as 




FIG. 1: Honeycomb lattice, graphene tt and vr* bands, and dispersion relation close to the K± points of the first Brillouin zone. 



with /(k) = 7o e*'^ '''" (e„i are the vectors that connect A-types site to the three nearest neighbor B-type sites, 
see Figure [H left panel) and 70 ~ —2.7 eV. It is readily shown that 

E^{K,ky) = ±70 [3 + 4cos(^^) cos(^) + 2 cos(fc,a)] ^'^ (2) 
(see Figure [J, central panel). 

The underlying symmetry and sublattice degeneracy of the honeycomb lattice can be unveiled from a different 
starting point. In a tight-binding representation, one can separate the two sublattices as 

3 

^ = E E ^R+e„-4R + /l.C, (3) 

R G B.L.m = 1 

by introducing creation {J^ and annihilation (-4,S) operators for electrons on site A or B, respectively. To 
diagonalize the Hamiltonian, one applies Bloch theorem to the Bravais lattice, which yields 

= 7^ E ^^""-^Ic ^Ue„ - 7^ E (4) 

with Nm the number of unit cells. By means of this transformation 

■n^ ^kBlAu, (5) 

kelBZ 

with 

3 
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In the sublattice basis, at a given crystal momentum, this also reads 



k£lBZ 





hk \ 


(4) 


{hk 


) 





Now, by performing a linearization close to K points, i.e. k = K± + p/fi. at first order in p — iftV, one finds 



'Hk+ = VF 



Px + iPy 
- ipy 



(7) 



(8) 



where vp ~ \/3a|7o|/2fi. is the Fermi velocity. Therefore, the effective Hamiltonian takes the form of a Dirac Hamil- 
tonian for massless particles 



(9) 



with Cri-.y^z the Pauli matrices. Finally, by identifying each of the K-points with the index, these Hamiltonians can 
be rotated into a diagonal form with the unitary operator {9 — arctan(pj,/pa;)) 



so that 



^lip)'HaipPi{p) ^ avF ( Q^' ]=-avF\P\crz 



(10) 



(11) 



The dispersion relation describes a cone- type structure in the reciprocal space (see Figure [TJ right panel), whereas the 
corresponding eigenvectors define a spinor (pseudospin) Q 



_1_ ( se'«/2 
V2 I e-''/2 



(12) 



The spinor is also an eigenstate of the helicity operator e = (l/2)(T.p/|p|, such that e|^'+(s = ±1)) = ±(l/2)|^'+(s)). 
Thus, eigenstates have a well-defined helicity (good quantum number). Note, however, that the above derivation is 
vaHd in first order in p, whereas trigonal warping, which appears at second order will bring substantial modifica- 
tions of band structure with consequences on transport properties. 

The electronic properties of the 2D graphene can thus be described by an effective massless Dirac fermion model in 
the vicinity of the CNP, with linear dispersion and electron-hole symmetry. These properties derived close to K point 
are also valid for ID systems such as metaUic nanotubes and wide armchair nanoribbons. However, other symmetries 
result in semiconducting systems with varying gaps. Semiconducting nanotubes and ribbons with increasing diameter 
(or width) show a linear downscaling of their associated energy gaps. By using proper boundary conditions, the 
electronic band structure of both types of systems can be analytically derived. 

In the presence of a uniform magnetic field B threading the 2D graphene sheet along the orthogonal z direction, we 
make use of the minimal substitution in the Hamiltonian. If we choose the first Landau gauge, the vector potential 
is A = {By, 0, 0) and the new momentum operator is 



Px 



-By 



(13) 



Py 



and so 



Let us define the following operators 



V2 h 



{tTx + ilTy) 



with 



J_io 

V2 h 



(14) 



(15) 
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They can be identified as annihilation and creation operators, since [77,77^] = 1. The Hamiltonian around K+ is now 

CO ?/ \ t,y 

with /3 = V2 — ^ . (16) 

Let us consider the square of the Hamiltonian 

nl^^P' U/3M . (17) 

V Tyyyt / \ i^^tj + 1 / 

We can identify rj^r] as the number operator, therefore the eigenvalues ofH (given by the square root of the eigenvalues 
of H^) are (H 

En =±PV^ ~ 31.65 VBy/7i with n = 0,1,2,3... , (18) 

where the eigenvalues E„ are expressed in meV and the magnetic field B in T. In contrast to ordinary two-dimensional 
electron gases, in graphene Landau levels are proportional to the square root on the magnetic field and to the square 
root of the integer number n. Moreover, a double degenerate zero energy level is present. 



B. Graphene nanoribbons band structure: nearest-neighbor tight-binding model 

Graphene nanonoribbons are strips of graphene that can be obtained by cutting a graphene sheet along a certain 
direction. Depending on this direction, the edges of the ribbon can be armchair-like (Figl2lla)) or zigzag-like (Fig. 
[Sl^a)). In the fohowing, we wih refer to an armchair ribbon composed of N dimers lines as iV-aGNR, and to a zigzag 
ribbon composed of N zigzag fines as 7V-zGNR. 

The band structures of ideal GNRs with width below ~ 100 nm and well defined ed g e s y mmetries (zigza.g or armchair 
types) are dominated by confinement effects and van Hove singularities (26l . [27l . [IJ [la, 113, Ell, [32, [33, [31 , similarly 
to carbon nanotubes. As for 2D graphene and carbon nanotubes, a nearest neighbor tight-binding model is found to 
describe low energy properties with a degree of approximation high enough for manyapplications. This a ppr oach has 
been widely employed for studying transport properties in pure or defected GNRs |26l . [27l . [H, [s^, [H, [33. l35l. [s^ . [37| . 



1. Armchair nanoribbons 



An armchair GNR unit cell contains A'' A-type atoms and TV B-type atoms, see Fig. [2lja). Thanks to the bipartite 
lattice of graphene, the total wavefunction of the system can be constructed by a hnear combination of A-type '^a 
and B-type ^'s sublattice wavefunctions. By applying Dirichlet boundary conditions for the wavefunctions at both 
edges of the ribbon {qy = 2pTr/{V3a{N + 1)), _p = 1, 2, .., A^), one gets [H 

N N 

(1/AA)|*) =. caY.Y.^''^^^' sin(V3(Zj,aj/2)|^,) +CB^^e'''="^^ sui{Viqya3 /2)\Bj) , (19) 

j = l xAj j=l xBj 

with N = yj2/Nx{N + 1) a normalization factor. Then, by rewriting the Schrodinger equation as a 2 x 2 matrix, 
eigenvalues and eigenvectors are readily found: 

S(fc,gy) = ±|7o(2e*'=''/2cos(%/3(?W2) + e"''^'')| =±|r(fc,gj,)| (20) 

and 

I*) = l/V2(^A ± y^r*(fc, qy)/T{k, qy)-^B^ ■ (21) 

By analyzing r{k,qy), we can see that the ribbon is metahic if = 3m + 2, where m is an integer number, and 
semiconducting in the other cases. In particular, we can obtain the value evaluate the energy gap Ajv as a function 
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of the N dimer chains: 



^3m 



ho\ ( 4 COS 

I70I ( 2 - 4cos 




3m + 1 

7r(m + 1) 
3m + 2 



(22) 



with Asm > A3„i-|_i > A3„i+2 = 0. However, as reported in section III C\ ah initio calculations do not predict the 
existence of metallic ribbons and a give a different gap size hierarchy, as discussed in Sect III Cl 

Further insight into the electronic structure of graphene nanoribbons can be gained by representing the Hamiltonian 
on the Bloch sums. If we order the basis as Ai, ^2, A3, Ajv_i, i?jv and -Bi, ^42, -B3, i?Ar_i, ^at, the Hamiltonian 




FIG. 2: (a) Structure and primitive cell of an armchair ribbon composed of A*' dimer lines, (b) Equivalent coupled chains 
with N orbitals each, (c) Energy bands for a metallic ribbon with A*' = 53. (d) 2D projection (at the F point) of the two 
wavefunctions just below and above the Fermi level for a 20-aGNR. (e) The same for a 35-aGNR. Blue or red colors correspond 
to opposite sign of the wavefunction. 
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can be split into four N x N blocks 



/ 1 



d ... \ 



1 1 



d* 



1 



d 



H(fc) = 70 



d* 







1 



(23) 



d 



1 1 



d* 



1 



V 



where d{k) ~ g^fc^/a rpj^^ diagonal blocks are tridiagonal matrices with diagonal elements and 70 off-diagonal 
elements, the off-diagonal blocks are diagonal matrices with alternating 70 cxp{±ikb/2) elements. Hamiltonian (|23|) is 
equivalent to the Hamiltonian of two coupled chains with N orbitals, as indicated in Fig. EJJb). Within each chain, the 
nearest neighbor orbitals are coupled by the hopping parameter 70 and have vanishing onsite energies. The interchain 
hopping is given by the parameters 70 exp(±zfc6/2). By diagonalizing l(23|) . the energy bandstucture is obtained. As 
an example, the energy bands of a metaUic armchair ribbon with = 53 are reported in Fig. [2b- Note the typical 
Dirac-like linear dispersion around k = 0. The mapping of the six Dirac points of 2D graphene into fc = can be 
easily understood by projecting the graphene band structure onto the axis corresponding to the armchair orientation 
[3^ 1. At fc = 0, the Hamiltonian (|23|) is real and our coupled chains turn into a ladder with tw o leg s and N rungs and 
all hopping parameters equal to 70. The eigenvalues of such a system are analytically known [39l. l40| and give again 
the energy gap An reported in Eq. l(22|) . 

The n-n tight-binding bandstructure of armchair nanoribbons with TV = 3m, 3m + 1 and 3m + 2 is summarized in 
Fig. Hfor N =9, 10 and 11. 

From the viewpoint of spatial symmetry, aGNRs, and GNRs in general, are very different from nanotubes. In fact, 
an ideal (perfect) carbon nanotube presents a well defined symmetry with respect to a large number of mirror planes 
containing the tube axis. As a consequence, the parity (even or odd) of the wavefunction has been demonstrated 
to have an impact on the transport properties of defected carbon nanotubes-based systems, as scattering can only 
occur between eigenstates with the same parity [ill . [H, 0, \4^. On the contrary, wavefunctions in GNRs do not 
always present a well defined parity associated to mirror refiections with respect to the axis of the ribbon. The spatial 
symmetry depends on both the edge termination (zigzag or armchair) and the even or odd number of chains that 
compose the ribbons. In the case for odd-indexed armchair ribbons (as the 9-aGNR and 11-aGNR in Fig. |4]), the 
system is invariant under mirror refiection with respect to a plane perpendicular to the ribbon and containing its axis. 
In the case of even-indexed armchair ribbons (as the 10-aGNR in Fig. |4|, the system is "asymmetric," in the sense 
that it is no more invariant under the mirror reflection operation. Still, it is invariant under glide plane symmetry, 
which consists of the mirror reflection followed by a fractional translation along the ribbon axis. Figures [2ld,e) show 
the 2D-projection at the F point of the two wavefunctions just below and above the Fermi level for perfect asymmetric 
armchair ribbons (calculated with the SIESTA package [45|). As clearly visible in Figs. [2l^d,e), the existence of a 
mirror symmetry plane in ribbons leads to a well deflned parity of the wavefunctions with respect to that mirror plane. 
The well deflned (or not) parity of the wavefunctions does not have any effect on the electronic transport properties 
of ideal unperturbed ribbons, but it may have a huge impact in the presence of defects. For example, the electronic 
transport properties of boron doped armchair ribbons have been proved, by means of ab initio calculations, to strongly 
depend on the symmetry of the ribbon, as B-induced potentials that preserve the parity of the wavefunctions do not 
affect the conductance of odd-indexed ribbons at low energies [Biel, B.; Blase, X.; Triozon, F.; Roche, S. submitted]. 
This suggests that scattering by certain defects might be suppressed, provided that the defects preserve the underlying 
symmetric geometry of the ribbon. 



For zigzag GNRs, the handling of the boundary conditions is slightly more involved. A possible approach is solving 
the Dirac equation requiring that the wavefunction vanishes at the edges of the ribbon [3l|. As evident from Fig. 
[3ja), the lower edge is entirely composed of A-type carbon atoms, while the upper edge is made up of B-type carbon 
atoms. Therefore, the boundary conditions can be imposed on the two sublattices separately. In particular, the 
wavefunction of the A-type sublattice is required to vanish on the upper edge and the wavefunction of the B-type 



2. Zigzag nanoribbons 
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sublattice is required to vanish on the lower edge. With these restrictions, the eigenfunctions of the Dirac equation 
can be separated in two groups. The first group includes states with wavenumber k > 1/W, where W is the width of 
the ribbon. They are surface states and decay exponentially from the edges of the ribbon as exp{—kx). The second 
group of eigenfunctions with k < 1/W corresponds to confined states with nodes along the transverse section of the 
ribbon. To obtain the exact shape of these bands within the simple tight-binding approximation, a different approach 
can be used, viewing the GNR as a macromolecule with a basis of four carbon atoms [i^ . 



(a) (d) 




FIG. 3: (a) Structure and primitive cell of a zigzag ribbon composed of A'" carbon lines, (b) Equivalent finite chain with 2N 
orbitals. (c) Energy bands for = 32. (d) 2D projection (at the F point) of the two wavefunctions just below and above the 
Fermi level for a 9-zGNR. (e) The same for the 10-zGNR. Blue or red colors correspond to opposite sign of the wavefunction. 

In this case too, further information can be obtained by expressing the tight-binding Hamiltonian on the Bloch 
sums. It is easy to see that the /c-dependent Hamiltonian is a tridiagonal 2N x 2N matrix with alternating off-diagonal 
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elements |42 



H(k) = 70 



/ 2c(fc) 
2c{k) 1 

1 2c{k) 

2c(/c) 



\ 





V 







... 2c(fc) 
2c(fc) / 



(24) 



where c{k) = cos(fca/2) and the Bloch sums have been ordered as Ai, Bi, A2, An , Bn ■ This Hamiltonian is 
equivalent to the Hamiltonian of the 2iV-orbitals chain reported in Fig. [3ljb). By diagonalizing matrix l(24|) numerically, 
the band structure is obtained for N = 32, see Fig. [3ljc). The energy bands present some typical features: (i) the six 
Dirac points of the 2D graphene are mapped into k = ±27r/ (3a), (ii) there are two partially fiat degenerate bands with 
energy between the Dirac points and the border of the Brillouin zone, the corresponding states are mainly located 
at the edges, (iii) the bands are highly degenerate at the borders of the Brillouin zone, (iv) at the Dirac points and 
close to the charge neutrality point the levels are equispaced. All these characteristics can be understood by simple 
considerations on the equivalent chain of Fig. \3lh). For example, at the edges of the first Brillouin zone we have 
c{k = ±7r/o) = 0, thus the 2A^-orbitals chain is split into two isolated orbitals at its end and iV — 1 couples of dimers. 
This configuration generates two eigenstates with zero energy and located exactly on the edges of the ribbon, and 
— 1-fold degenerate eigenvalues E = ±70, as observed at point (iii). If we consider a very large ribbon, our chain 
can be thought as semi-infinite. In this case, an analytical expression for the retarded Green's function projected on 
the first (Ai) site can be obtained ji^l 



G{E,k) = 



Ajy + 7g ± ^{E^ - 47gc2 + ~ ^l',E^ 
HE 



(25) 



where the sign in front of the square root is positive if < c(fc) < 1/2 and negative if 1/2 < c(fc) < 1. In order to 
have a nonvanishing density-of-states (DoS), the discriminant under the square root must be negative, since the DoS 
is proportional to the imaginary part of the Green's function. Therefore, the four curves E(k) = ±70 ± 270 cos(fca/2), 
along which the discriminant is vanishing, define the boundaries of the band structure and their intersection at 
k ~ ±27r/(3a) identifies the folded Dirac points. As in the case of armchair nanoribbons, the mapping of the Dirac 
points could be predicted by the projection of the 2D graphene bands onto the axis corresponding to the zigzag 
orientation. When E ^ Q, the retarded Green's function has a continuum of poles for < c(fc) < 1/2, which 
correspond to the fiat bands for |fc| > 27r/(3a). The structure of the energy bands around the Dirac points and for 
finite width ribbons is obtained by considering that cik = ±27r/3a) = 1/2, i.e. the hopping energies of the equivalent 
chain become all equal to 70. The eigenvalues of such an Hamiltonian are known in the literature [i^ and, for large 
ribbons and around the CNF, can be approximated as 



En = I70I 



A^ + 1/2 



with 



0,±1,±2,. 



(26) 



As observed in point (iv) , this means that these levels are equispaced and the spacing is inversely proportional to the 
width of the ribbon itself. 

The n-n tight-binding bandstructure of zigzag nanoribbons are summarized in Fig. |4]for A^ =10. 

As in the case of aGNRs, the even or odd number of chains that compose the ribbons determines the mirror or 
glide plane symmetry of zGNRs. Even-indexed zigzag ribbons (as the 10-zGNR in Fig. |4]) are invariant under mirror 
reflection with respect to a plane perpendicular to the ribbon and containing its axis. Odd-indexed zigzag ribbons 
are invariant under glide plane symmetry. The parity of the wavefunctions with respect to the spatial symmetry of 
the GNR is clearly visible in Figs. [3td,e), where we show the 2D-projection at the T point of the two wavefunctions 
just below and above the Fermi level for perfect asymmetric (d) and symmetric (e) zigzag ribbons (calculated with 
the SIESTA package [Hi). The spatial symmetry of zGNRs and the respective parity of the wavefunctions play an 
important role in determining the electronic transport properties in the presence of external fields (iol . Isol . Isil . Is^ l . For 
example, a valley- valve effect (blocking of the electrical current by a p-n junction) has been demonstrated theoretically 
in zigzag ribbons in the presence of an external potential j4§| . This effect depends critically on the parity of the number 
of zigzag chains across the ribbons [H3| and has been successfully explained by reaHzing that, in case of even-indexed 
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FIG. 4: Band structures of 2D graphene, (9,0) and (10,0) zigzag and (10,10) armchair CNTs, and 10-zGNR, 9-10-11 aGNRs 
obtained by n-n tight-binding model, third n-n tight-binding model and ab initio calculations. The energy bands of 2D graphene 
are plotted along the K-F-M-K' direction. The energy bands of the GNRs and CNTs are plotted in the one-dimensional first 
Brillouin zone. The wavevector k is in unit of tt over the lattice constant. 



zGNRs, the symmetry of the ribbons is preserved by the considered superimposed potential, thus forbidding potential- 
induced transitions between opposite parity eigenstates [131 ■ A similar explanation has been found for the case of 
an applied bias in zigza g ri bbons, when transitions between the valence and the conduction bands are only allowed 
for asymmetric ribbons [5l| , and the opening of a conductance gap in the vicinity of the Fermi level is expected for 
symmetric zigzag ribbons. This intriguing particularity has no analogue in nanotubes. 
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3. Carbon nanotuhe hand structure 

As for graphene nanoribbons, the electronic properties of nanotubes are strongly modulated by small structural 
variations, in particular, their metaUic or semiconducting character is determined by the diameter and hehcity (chi- 
rality) of the carbon atoms in the tube. This dependence is easily understood starting from the energy dispersion 
relations of graphene around the Fermi level, by retaining only the tt-tt* bands. Consequently, rolling the graphene 
into a cylinder imposes periodic boundary conditions for the wavefunctions along the circumference of the tube and 
results in a quantization of the momentum component along this direction. For derivation details of the tight-binding 
graphene dispersion relations, and the subsequent zone-folding procedure to obtain the band-structure of nanotubes, 
we refer the reader to Ref. The obtained bands for armchair tubes with helicity {N, N) are, for instance, given by 

(fc) = ±70 y 1 ± 4 cos ^ cos ^ + 4 cos^ ^ , (27) 

where q (= 1,2, ...,27V) specifies the discrete part of the wavevector perpendicular to the tube axis (i.e., the band 
index), k is the continuous component that describes eigenstates in a given sub-band (— tt < ka < tt) and a = 2.46 A 
is the graphene lattice constant. 

The bandstructure for the armchair (10,10) and the zigzag (9,0) and (10,0) nanotubes calculated within the n-n 
tight-binding model are shown in Fig. [H 

C. Ab initio results 

In spite of the good description the most simplified n-n 7r-orbital tight-binding model provides for the electronic 
properties of 2D-graphene and carbon nanotubes, some particular features regarding mainly the effect of edges in 
graphene nanoribbons have only been unveiled by means of ah initio simulations pSi . IH^ . Issl . IsgI . IstI . [58| or more 
sophisticated tight-binding models [H^, [13, [fiH. The stability of edge states in zigzag GNRs has been analyzed both 
with and without H atoms by means of ah initio simulations. The almost fiat bands in the vicinity of the Fermi 
level, originated by the presence of very localized states at edge atoms of the ribbon, become more dispersive, thus 
producing a peak in the conductance of the system, that jumps from one to three conduction channels in a region of 
w 0.3 eV below the CNP of the ideal ribbon. The dispersion at the edge state predicted by first-principles calculations 
has been attributed to the interaction with next-nearest neighbors, that decreases the energy eigenvalue of the edge 
state [gOI, though there are still some discrepancies with respect to experimental data [62 . £3,, i6J]. However, the 
metallic character of zigzag ribbons is preserved as far as the spin degree of freedom is not taken into account [55| . 

On the other hand, ab initio studies have demonstrated that there are no truly metaUic armchair graphene nanorib- 
bons ([H, [13, E^l). Even for those armchair ribbons predicted to be metaUic by the n-n tight-binding model, a gap 
opens, thus modifying their character from metallic to semiconducting. The magnitude of the gap, however, decreases 
with increasing ribbon width, and for a ribbon of w 5 A with a metaUic behavior predicted by the n-n tight-binding 
model, the ab initio estimated band gap is only « 0.05 eV. The magnitude of the gaps of the semiconducting ribbons 
predicted by the n-n tight-binding model also differs from those estimated by means of first-principles calculations. 
There have been several explanations for the origin of these gaps, supported by more sophisticated tight-binding 
models that take into account edge distortion [13,11§|, up to three nearest neighbors interactions [13|, or both [ll|. 

The ab initio bandstructure of 2D graphene and different types of nanoribbons and nanotubes are shown in the last 
column of Fig. |4l All the ab initio calculations have been performed with the SIESTA package [45| . 

D. Third nearest neighbors interaction tight-binding model 

1. Changes in 2D graphene band structure 

In section III Al it was mentioned that the band structure of 2D graphene can be calculated using the nearest- 
neighbor tight-binding approximation [l^. This leads to Eq. ([2]), which is plotted in Fig. [1] (central panel) and Fig. 
m We will now see how this result is modified by taking into account higher orders of interaction and orbital overlap. 

Figure [UJa) graphically shows the distances between first, second and third nearest neighbors. These distances 
are given by di = 1.42 A, 6,2 = V^di « 1.73di and = 2di. From Fig. [Uja), we also see that every atom in 
the graphene lattice has three nearest neighbors, six second-nearest neighbors and three third-nearest neighbors. 
The small difference in length between the distances of second and third nearest neighbors (^2 and d^) suggests the 
inclusion of both of these interactions. Figure [5ljb) allows us to evaluate which interactions exist between the unit 
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X 



FIG. 5: (a) An atom in the graphene lattice with its first, second and third neighbors, indicated by three dashed circles. Second 
neighbors belong to the same sublattice of the considered atom, first and third neighbors to opposite sublattice type, (b) A 
unit cell of the graphene lattice and its adjacent unit cells. Boundaries between unit cells are marked in green, they are labeled 
in blue. Purple lines represent first, orange lines for second and red lines for third nearest neighbor hopping. 



cell here labeled (0,0) and its neighboring unit cells. We can therefore write the tight-binding Hamiltonian, including 
up-to-third nearest neighbor interaction: 

H(fc) = Ho,o + e"=-''iHi,o + e-"=-"iH_i,o + e"=-"^Ho,i + e-"=-"^Ho,_i + e"=-(°i+°^)Hi,i + e-"=-^"i+°^)H_i,_i 

+ e'^^-f^i-^^^Hi^^i +e-''=-(''i-''^)H_i,i (28) 

where ai ~ (\/3, 3)a/2 and a2 = (V3,-3)a/2 are the translation vectors of the hexagonal lattice and 



\ --yo eg ; -71 / V -71 



-12] „ _ (-11 \ _ / 

) ^"-1 - 1-70 -71 j 1-72 



(29) 



with = for i 7^ j. First, second and third nearest neighbor hopping constants are here called 70,71 and 72, 
in accordance with Fig. [5ljb). The onsite energy is indicated by eo- Analogously, we define ctq, cti and (T2 to be the 
overlap matrix elements between first, second and third nearest neighbors, respectively, and we get a formula for S 
that looks just Hke Eq. l(28l). with every H replaced by an S and 



G _ ( 1 cro\ o ^ f^i '^2\ Q _ /^fl CTO 

*o,o — \ „ 1 I *i,-i ~ \ rr /X ) ''o.i 



« /O <J2\ „ _ fai 0\ f 

^1,1 - n n =>0,-l - I _ „ ^-1,-1 - 



(30) 

,0 ; ^-'-^ - {ao a,j - \a2 ()] 

where of course Si^- = Sj^i for i ^ j. We now define 
go(fc) = l + e'^°-^ +eifca2 

gi(k) = e""'"! + g-ifc-oi _|„ gifc-a.2 _|„ g-ifc-a2 _j_ gifc-(oi-02) _|_ g-ifc-(oi-02) {^^) 
g2[k) = gi'=-("i+'*2) _|_ gifc- (0.1-12 ) _|_ g-ifc-(ai -12) ^ 

which permits us to write the Hamiltonian H as 

TT_/' eo-7i.9i(fc) -7o5o(fc) - 72.g2(fc)A /o9^ 
I -7050* (fc)- 72.92 (fc) eo-7i5i(fc) j ^ ^ 



and the overlap matrix S as 



c_( 1 + 0-151 (fc) cro.go(fe) + 0-252 (fc)\ /oq^ 
* - Uo5o(fc) + ^2.92* (fc) 1 + ^151 (fc) ' ■ ^ ^ 
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Now all that is left to do is to solve the generalized characteristic equation 

|H-£:s| = 



Haa — ESaa Hab — ESab 
^AB ~ ^^AB EIaa — ESaa 



= {Haa - ESaa? - {Hab - ESab){H*ab - ES\b) = 

(34) 



^ E-^iSl^ - l-^AsH + E{HabS*ab + H*abSab ~2 HaaSaa) + H^^ - \Hab\^ = , 
where we introduced ei.2,3,4 as a shorthand notation. From comparison between Eqs. (|32l) and l(33|) we see that 

ei = (1 + c^i,9i(fe))^ - |cro.go(fc) + 0-252(^)1^ 
£2 = -[(7o5o(fc)+7252(fe)) X (cro5(5(fe) + cr2.g2('s)) + c-c.] 



(35) 

£3 (eo - 7igi(fc)) X (1 + cri.gi(fc)) 

£4 = (eo - 'yigi{k)f - |7o5o(fc) + 72.g2(fc)|^ . 



We are thus left with the following dispersion relation for 2D graphene: 



^± ^ -(£2-2£3)±v/(e2-2£3)^-4£i£4 ^3g^ 

2£l 

The values for 70-2 and (T0-2 have been calculated by Reich et al. by comparing and fitting band structures from 
ab initio and tight-binding calculations [qBI- The third-nearest-neighbor tight-bindin g b and structure with these 
parameters is plotted in Fig. IH Recently, more sophisticated GW ab initio calculations [63| for bilayer and few-layer 
graphene provided ~ 20% larger values of 7s. This is in better agreement with some experimental measurements, in 
particular angle-resolved photoemission (ARPES), where the long-range Coulomb interaction plays an important role 
and the band structure at high energies requires an more accurate description. 



2. Effects on the band structure of GNRs 

It was already mentioned that the metallicity of armchair GNRs of widths = 3m + 2 (m e N) is only found in 
the n-n tight-binding approximation, while ab initio calculations show all armchair GNRs to be semiconducting (cf. 
section [TTB] and section ITTC|) . This finding is reproduced by third-nearest-neighbor tight-binding calculations. The 
gap opening can be quantified by comparing n-n and 3n-n tight-binding results with each other (cf. Fig. [6]). The 
relative opening or clo sing of the ban d g, ap decreases with the width as 1/N . Furthermore, it has been shown that 
edge deformations (33l . l34| and strain [63, |6^ significantly change the band gap of armchair GNRs, turning metallic 
ribbons into semiconducting ones, even when calculated with the n-n tight-binding model. The metallicity of certain 
armchair ribbons should therefore be considered an artifact of the n-n tight-binding approximation and not a robust 
feature of these ribbons. 

In the case of zigzag GNRs the 3n-n tight-binding model shows robust metallicity and a higher dispersion of the 
edge state than in the case of only nearest-neighbor interaction (cf. Fig. |4]) . The full ab initio results show a band gap 
that opens due to magnetic ordering (H^]. This feature is not reproduced by the tight-binding model employed here, 
but one can account for such spin- related effects by including a Hubbard term in the Hamiltonian [s^ . ItoI . ItiI . It^ I . 

It can thus be said that the third-nearest-neighbor tight-binding method leads to significantly better results than the 
nearest-neighbor method when compared to ab initio results, while keeping within the simpHcity of the tight-binding 
model. 

The 3n-n energy bands for 2D graphene and different types of nanoribbons and nanotubes are summarized in the 
fourth column of Fig. [4l 



E. Transport in nanotubes and nanoribbons 



When neglecting disorder effects, a carbon nanotube of length L in between metallic contact reservoirs is a ballistic 
conductor with a L-independent conductance given by the energy-dependent number of available quantum channels 
N±{E). Each channel carries one conductance quantum Gq = 2e^//i, thus G{E) = 2e^//i x N±{E), including spin 
degeneracy. This only occurs in case of perfect (refiection-less) or ohmic contacts between the CNT and metaUic 
voltage probes (73l. [74|. In this regime, the expected energy-dependent conductance spectrum is easily deduced from 
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FIG. 6: (a) The band gaps of armchair GNRs as a function of their width A*" in n-n and 3n-n tight-binding. Because of the three- 
fold periodicity in the band gap's width, it is convenient to plot three different curves or sets of points for A*" mod 3 = 0, 1, 2. 
The symbols show n-n, the lines corresponding 3n-n results, (b) The gap opening/closing obtained from a comparison of n-n 
and 3n-n tight-binding results. The symbols show the opening (or closing) of the band gap in armchair GNRs of the respective 
widths. The lines are fit to a/N, where a is the fitting parameter. This is analog to the case of zigzag carbon nanotubes in 
which the gap decreases as l/d, where d is the diameter of the tube. 



band structure calculations, by counting the number of channels at a given energy. For instance, armchair nanotubes 
present two quantum channels at the charge neutrality point, which result in G{Ef) — SGq. At higher energies, the 
ballistic conductance increases as more channels become available to conduction. 

In contrast, the eigenstates of defect-free armchair GNRs in the first plateau are constrained by the hard- wall 
(Dirichlet) boundary conditions which require a node at the edges, thus excluding the cosine solutions [IHl- As a 
result, a singly degenerate band is left in the first conductance plateau. 

Within the linear regime near the CNP, a closed expression for the the number of channels is given by 

7V,(i.).2 + 4(^|£;-£;H) (37) 
in the case of an armchair CNT. For CNTs of other chiralities and GNRs, similar expressions can be obtained. 



1. "Pseudodiffusive" transport in clean graphene nanoribbons 

Transport through short and wide clean graphene monolayers around the charge neutrality point shows very peculiar 
behaviors analogous to those observed in diffusive coherent disordered conductors. This turns out to be a mere but 
extremely surprising coincidence, due to the particular distribution of the transmission coefficients. In two first 
theoretical studies, Tworzydlo et al. [iHl and Katsnelson [t^ have investigated the transmission coefficients for a 
ribbon of length L and width W in the fimit of short and wide armchair ribbons (L < W) with hardwall or smooth 
confining potential at the edges. The system is kept at energies around the neutrality point, and it is connected to 
two leads at high potential with a large number N of active conductive channels. The electronic transmission is thus 
mainly sustained by tunneling through evanescent modes, and an analysis based on the Dirac equation leads to 



(38) 
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for the hardwall confining potential, and a very similar expression for the smooth potential [75||. As a consequence, 
in the limit of a very large number of transmission channels and high aspect ratio {N » W/L oo) 

where G is the conductance and F is the Fano factor, i.e. the ratio between shot noise and Poissonian noise. In this 
limit, the system is ohmic with minimum conductivity fTmin = 4e^/7r/i and Fano factor equal to 1/3, as in the case 
of a diffusive conductor. This theoretical result has been tested numerically by calculations based on a tight-binding 
model (tsI. IttI. and a good agreement was found, see Fig. [71 
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FIG. 7: Conductivity cr of a zigzag ribbon as a function of its aspect ratio (length L over width W), in the case of graphene 
leads and energy just above the charge neutrality point. For short and wide nanoribbons (0.1 < LjW < 0.4), the value of the 
conductivity is very close to Ae^ /-rrh. This figure is taken from [Cresti, A.; Grosso, G.; Pastori Parravicini, G. Phys. Rev. B 
2007, 76, 205433]. 

Schomerus has recently proved [tH] that this behavior is essentially independent of the specific configuration of the 
contacts, provided that the graphene system is very close to the neutrality point and the leads support a large number 
of propagation modes. The introduction of disorder is expected to affect the Fano factor, but the literature does not 
agree on the results. In one case [t^ the Fano factor turns out to be lowered by disorder, but still keeps a value close 
to 0.3. In another case [80|, the Fano factor is found to be a little lower than 1/3 for very weak disorder with a small 
peak at the neutrality point, and then to stay around 1/3 for short and wide ribbons for weak disorder, and to go 
over 1/3 for stronger disorder. 

Three recent experimental works [sU - ls^ -lssI on exfoliated graphene samples confirm the main theoretical predictions. 
Miao and co-workers j8l| have first evidenced that for wide and short graphene strips {W/L > 4) the minimum 
measured conductivity approaches the value ie'^/irh. This result has been observed in particular for very broad 
ribbons {W >>IA, because in this case boundary effects are negligible. DiCarlo and co-workers [s^l and Danneau 
and co-workers [83| performed measurements of the shot noise and observed the predicted value = 1/3 for clean 
samples. Moreover, some other predictions [80| seems to be confirmed, i.e. the increasing value of the Fano factor 
for weakly disordered systems [83| and the presence of a (larger than expected) peak at the neutrality point for very 
clean ribbons [83 |. 

The coincidence between the behaviors of coherent disordered conductors and short and wide ballistic graphene 
ribbons is a puzzHng issue, whose origin, if any, is not clearly understood. 



2. High magnetic fields and spatial chirality of currents in GNRs 



As seen in section [TTAI a homogeneous perpendicular magnetic field B that threads a 2D graphene sheet induces 
Landau levels spaced as ^/n and proportional to \fB (25l . [83 | . A prescription for the inclusion of magnetic field effects 
into the band structure of carbon-based nano-networks can be found in Ref . [s^ . An energy structure similar to that 
of 2D graphene is also observed for bulk electrons in nanoribbons, provided that the width W of the ribbon is larger 
than the magnetic length. At the borders of the Brillouin zone, the bands bend upward or downward due to the finite 
width of the ribbon. The corresponding magnetic states are located at the upper and lower edges. 

In the simple n-n tight-binding model described in section III B[ we can account for the magnetic field by means 
of the Peierls phase factors on the hopping parameters. By choosing the gauge properly, the Hamiltonian preserves 
the translation invariance along the longitudinal axis of the ribbon and can be conveniently expressed on the Bloch 
sums basis. Armchair and zigzag ribbons are still equivalent to a 2A^-orbitals chain and two coupled N orbital chains. 
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FIG. 8: (a) Coupled chains equivalent to an armchair ribbon threaded by a perpendicular and homogeneous magnetic field B. 
(b) Energy bands of an aGNR with A'' = 611, corresponding to a width of about 75 nm, for B = 10 T. (c) Chain equivalent to 
a zigzag ribbon threaded by a perpendicular and homogeneous magnetic field B. (d) Energy bands of a zGNR with A'' = 353, 
corresponding to a width of about 75 nm, for B = 10 T. 



respectively, but the hopping energies now depend on the chain index n = 1,2, N [43|. In the case of aGNRs, the 
coupHng parameters between the two chains become (see Fig. EJa)) 



dn (k) = exp 



, kb na N +1 



(40) 



where a = 27r$(i?)/$o is proportional the ratio between the magnetic flux $ through a hexagonal plaquette of the 
honeycomb lattice and the elementary magnetic flux $0. In the case of zGNRs, the hopping fc-dependent parameters 
between the orbitals of the chain become (see Fig. [8l[c)) 



c„(fc) 



cos 



ka 
T 



na 



N + l 



(41) 



A direct numerical diagonalization of the Hamiltonian provides the energy bands for a iV = 611 armchair ribbon 
(Fig. [8l[b)) and N — 353 zigzag ribbon (Fig. [Hd)) for B = 10 T. We consider ribbons much wider (about 75 nm) 
than those considered in the absence of magnetic field, in order to have the onset of Landau levels at relatively small 
magnetic fields. 
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In the case of armchair ribbons, see Fig. [Hb), we can clearly observe the sequence of double degenerate Landau 
levels. At the borders of the Brillouin zone, the average transverse position of the states approaches the edges of the 
ribbon and the energy bands bend upward for electron-like particles and downward for hole-Hke particles. In fact, the 
particles cannot complete the cyclotron orbit due to the edges and as a consequence their energy rises. As in ordinary 
two-dimensional electron gas, the states on the right side of the Brillouin zone are located on the upper edge of the 
ribbon, while the states on the left side are located on the lower edge. The double degeneracy of the bands is removed 
for edge states. This can be explained in terms of the Dirac equation since the boundary conditions for armchair 
ribbons entail the admixing of the valleys [sHl- As a consequence, the corresponding wavefunctions hybridize thus 
giving rise to the observed structure. 

In the case of zigzag ribbons, see Fig. [HJd), we observe a positive and a negative set of Landau levels around each of 
the two valleys, and a zero energy level that extends from one valley to the other. The states on the right side of each 
valley are located on the upper edge of the ribbon, while the states on the left side of each valley are located on the 
lower edge. Again, when the average transverse position of the states approaches the edges of the ribbon, the energy 
bands bend upward for electron-like particles and downward for hole-like particles. The structure of the energy bands 
around the two valleys is exactly the same. Again, this can be understood by considering the Dirac equation. In 
contrast to armchair ribbons, the boundary conditions for the wavefunction do not mix the two valleys and therefore, 
close to the charge neutrality point, they behave independently of each other. The Dirac equation also explains the 
peculiar structure of the two lowest Landau levels. The corresponding wavefunctions are extremely localized at one 
edge for a sublattice and behave as ordinary magnetic states for the other. 

From the transport perspective, the first important consequence of the magnetic electronic structure of GNRs is 
the theoretical prediction and the experimental observation of the so called anomalous (of half integer) Hall effect 
0, [13, [13, [Hi, i.e. the magnetoconductivity is given by 

G=^(n+i). (42) 

The large separation between the first Landau levels allows the observation of this phenomenon at room temperature. 

In addition to this, magnetic edge states have peculiar and interesting properties. Since the direction of the carrier 
group velocity depends on the sign of their energy band slope with respect to the wavevector k, electron-Hke particles 
can only flow from the left side of the ribbon to the right side along the upper edge, and from the right side to the 
left side along the lower edge. In the case of hole-like particles the energy bands are reversed and they can only flow 
from the left side of the ribbon to the right side along the lower edge, and from the right side to the left side along 
the upper edge. The spatial separation between the conductive channels that support current in opposite directions 
is called spatial chirality. In the case of graphene nanoribbons, electron- and hole-like particles have opposite spatial 
chiralities This property may be exploited to manipulate the current and obtain a current switch [8§|. In fact, 
the edge the transport currents flow along is determined by the position of the electrons energy with respect to the 
charge neutrality point, and their energy can be varied by means of top or back gates, even locally. 

III. DISORDER EFFECTS 

Transport properties in graphene-based materials also turn out to be strongly affected by disorder, which can 
originate from impurities (charges trapped in the oxide, chemical impurities, etc.), topological defects (vacancies, 
edge disorder,...), or long range deformation modes (ripples) in 2D graphene. Importantly, the dominant scattering 
processes and resulting transport features are very dependent on the range of the disorder potential and the robustness 
or destruction of the underlying sublattice symmetries. 

For massless Dirac fermions, a long range scattering potential, i.e. with Fourier components ^(q) such that q ^ K, 
will strongly reduced the intervalley scattering probability between the two non-equivalent Dirac nodes (K_ K+). 
In the one-dimensional case provided by armchair CNTs, this results in a full suppression of backscattering probability 
as demonstrated by Ando and co-workers (90l . [9l] |. The impact of long range disorder in two-dimensional graphene 
remains more debated, with opened issues concerning how weak or strong localization regimes are genuinely affected 
by speciflc Dirac fermions properties [92, 93. 94.i95|. 

In contrast, the presence of short range disorder (q ~ K) will allow for all possibilities of intravalley and inter- 
valley scattering events between K_|_ and K_, leading to stronger backscattering and localization effects, although a 
true Anderson localization in two-dimensional graphene remains flercely debated theoretically, because of the "Dirac 
nature" of low-energy excitations [sB, [13, [l^] . 

Disorder effects in the quantum coherent regime can yield locaHzation regimes [13, llOOl . llOll . Il02l |. and indeed 
several experimental evidences for weak and strong Anderson localization regimes have been reported in disordered 
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carbon nanotubes 0, [Tol, [loi, [io3, \lM, [llil, [HI, but one should also stress that the effects of 

disorder on quantum transport in graphene-based materials of lower dimensionality, such as carbon nanotubes and 
graphene nanoribbons, are expected to be maximized if compared to the case of 2-dimensional graphene. Indeed, low 
dimensionality and confinement effects yield strong modifications of band structures with the appearance of van Hove 
singularities (vHs), close to which energy dispersion -or wavepacket velocity- is very low. In the forthcoming sections, 
we will focus on short-range disorder effects by using the Anderson-type disorder potential, and we will review its 
effects on electronics states, mean free paths and locaHzation phenomena. 



A. Density of states in weak disordered CNTs and GNRs 



The spectrum of a system may be affected by weak disorder only through small energy shifts in the energy. A 
constant density of states will therefore remain unaffected by weak disorder. Regions of slowly changing DOS will 
only show small effects. Discontinuities and van Hove si ngularitie s will, in contrast, be strongly affected and appear 
smeared out at the energy scale of the disorder strength [lid . Ill4| . 

The DoS of a general quantum wire under the influence of Anderson disorder can be obtained via an algorithm based 
on diagrammatic perturbation theory that takes into account localization effects by inclu ding multiple scattering [ll3l |. 
By dropping the crossing diagrams within the noncrossing approximation (NCA) |ll5l |. we can write the self-energy 
S (E) to all orders as a recursive expression, which can then be iterated numerically until self-consistency is reached. 
Though the a pplic ability of the NCA is not obvious, it can be justified by comparing the contribution of various terms 
at low orders |ll3| . 

The original formulation of this approach is restricted to the special case of CNTs where all atoms are equivalent 
through symmetry and the self-energy takes the same value for all atoms. It can, however, be generaHzed to arbitrarily 
structured quantum wires using matrix notation. The self-energy S (E) that accounts for disorder is a diagonal matrix 
obeying the recursive relation 



{E + iO+ -Ha -^{E))' 



(43) 



For a periodic system, the self-energy has the same periodicity as the Hamiltonian. The block-tridiagonal matrix 
{E + 10+ — Tip — S ( E)) can therefore be inverted numerically using a highly convergent renormalization-decimation 
algorithm [1161 Ill7| . allowing us to go beyond the energy range near the Fermi energy, where the special band 
structure allows analytic inversion. 

Starting with E = 0, each numerical iteration of this recursive relation is equivalent to one additional perturbative 
order. Typically, convergence is achieved after less than ten iterations, except for energies near a van Hove singularity, 
where hundreds of iterations may be necessary. This clearly indicates that low-order perturbation theory breaks down 
near band edges. 

The LDoS of each orbital i in the unit cell can now be obtained directly from the imaginary part of the Green 
function 



g,{E) = --Im 



{E + iO+ - Ho - {E)) ' 



(44) 



Figure [9] shows this quantity in direct comparison with the numerically exact value obtaine d by sample averaging. 
The slight deviation visible at the fianks of the van Hove singularities is caused by the NCA |113| . The elastic mean 
free path iei based on the DoS of a disordered system, displayed in Fig.[TOl is no longer a purely perturbative quantity, 
but it takes into account the scattering into localized states present at any given energy. 
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FIG. 9: Density of states g{E) of an infinite armchair (5,5) CNT under the infiuence of Anderson disorder. Dashed line: g (E) 
of the clean system displaying the van Hove singularities. Solid line: g {E) from Eq. I|44p . Data with error bars: values obtained 
numerically by sample-averaging. Adapted from [Nemec, N., Richter, K.; Cuniberti, G. New J. Phys. 2008, 10, 065014]. 



The DoS of aGNRs, see Fig. Elja), is obtained straightforwardly and shows the expected smoothing of van Hove 
singularities (including those at the gap edges) in complete analogy to the CNT. For zigzag-edge GNRs (Fig. [T4lfb)). 



19 




0.0 0.5 1.0 1.5 2.0 2.5 

^(eV) 

FIG. 10: The elastic mean free path lc.i in a (5,5) CNT with Anderson-disorder of strength cr^ = 0.5 eV. Both lines are 
obtained from Eq. I|49p . In the case of the dashed line, g{E) and A'ch correspond to the values of the clean system, resulting in 
discontinuities at the band edges, similar to Fig. ifTTjl . For the solid line, the disorder effects on g{E) and A^ch are taken into 
account using Eqs. I|43p and l|44p . 

the most prominent feature in the DoS is the strongly broadened edge state around the Fermi energy, whose shape 
displays a clear deviation of the NCA from the numerical exact calculation. 

B. Elastic mean free paths 

The elastic mean free path {lei) is a key quantity in mesoscopic transport. Hereafter we focus on the case of 
short range disorder, that allows us to illustrate common properties of transport length scales in graphene-based 
low dimensional materials. The Anderson disorder model is the most generic case for investigating locaHzation 
phenomena in low dimension. In this model, the onsite energies of Pz orbitals assume random values within an 
interval [—W/2,W/2] with a given probability distribution. Hereafter we assume a uniform probability distribution, 
i.e. V = l/W. 

In a situation of weak disorder, within the Born approximation scheme, iei can be derived at a certain degree of 
approximation, thus enabling the possibility to extract an analytical expression. The simplest approximation for 2D 
graphene can be derived as follows. The total density of states can be written as 

p{E) = (^a'/2n)\E\l{hvFf. (45) 

As a result, writing iei = vpT, and using a simple Fermi Golden Rule approach for the elastic scattering time r 
(r-i = {2-K/h)p{EF)W'^/ 12), we obtain 

iei - (7o/M^)'a|7o|/|^|, (46) 

which diverges when \E\ ^0. This crude estimation pinpoints a difficulty in calculating transport length scales when 
the Fermi level Hes close to the Dirac point. A numerical calculation within the Kubo approach allows the evaluation 
of ^ei at a quantitative level in 2D disordered graphene with Anderson scattering potential. In Fig. [19] (inset) (adapted 
from [ll8l |). £ei is shown for W = 1, 1.5, 2, 2.5 in 7o/2 unit and ranges from several tens of nanometers down to a few 
nanometers close to the Dirac point (for « 3.4 eV). 

In quasi-lD systems (such as CNTs and GNRs), scattering angles are restricted to two cases: Forward scattering 
events at zero angle lead to momentum relaxation but do not affect the elastic transport length scale. Only the 
backscattering events at an angle of tt are taken into account for the derivation of the elastic mean free path £ei ■ By 
using the Anderson disorder model. White and Todorov derived an analytical formula for the elastic mean free path 
{iei) close to charge neutrality point [9l|. For armchair metallic nanotubes, they obtained 

4i = 18V3acc(7o/W^)'A^, (47) 

showing that, for a fixed disorder strength, £ei will upscale linearly with the nanotube diameter, a property unique 
to these systems and pinpointing long baUistic systems. Triozon and co-workers [ll9l | numerically confirmed such 
prediction and further reported on the strong energy dependence of £ei close to the onsets of new subbands. Similarly, 
£ei was derived in metallic N-aGNRs (s^l as 

4i = 12(7o/M^)'(iV + l)acc. (48) 

Therefore, both low dimensional carbon systems show a mean free path that diverges with increasing diameter or 
ribbon width for a fixed disorder strength W . Note however that only armchair nanotubes really display ID massless 
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FIG. 11: The elastic mean free path lei in a (5,5) CNT with Anderson-disorder of strength = 0.1 eV. The numerical results 
are obtained from the transmission of a single disordered unit cell, averaged over 4000 samples as £ei — L/ (Nch/ (T) — 1). The 
analytical result is obtained via Eq. (|49|l using g{E) of a disorder-free system. 



Dirac fermions close to the charge neutrality point, since there is a gap opening in all GNRs due to edge boundary 
conditions. 

An expression for £ei that holds at arbitrary energies can be derived from the Fermi Golden Rule [ll4l | 



(49) 



with £uc the length of the unit cell and the sum running over all orbitals i within one unit cell. In this form, the 
expression can be appHed to arbitrary quantum wires, including GNRs, where it also covers the special case of edge 
disorder by making ct^. dependent on the orbital number i. 

Neglecting multiple scattering, the elastic mean free path £e and the diffusive transmission Tdiff are defined in 
terms of the LDoS gi {E) of the disorder-free system. Likewise, A^ch is defined by the leads, where it follows an exact 
integer step function. Near band edges, this diffusive transmission is discontinuous, as can be confirmed numerically 
to arbitrary precision, computing it as the sample average (T) of the transmission of many disorder configurations, 
as displayed in Fig. ifTTj) . 

For stronger disorder and extended disordered regions, the elastic mean free path is no longer a purely perturbative 
quantity due to the fact that it depends on the density of states that has to include non-perturbative effects near van 
Hove singularities. We can, however, retain Eqn. (|49|) by simply including the non-perturbative effects in gi{E) using 
Eqn. l(44|) . Furthermore, the number of channels A'ch must also take into account the non-perturbative effects near 
band-edges, which can be achieved by including the self-energy term Eqn . (|43|) in the calculation of the transmission 
of through a cross section of an infinite quantum wire as described in Ref. [ll4l |. Incorporating these effects, we obtain 
the elastic mean free path in an infinitely long disordered quantum wire as displayed in Fig. [TOl 



C. Quantum interference effects and localization phenomena in disordered graphene-based materials 

The knowledge of the mean free path £e in quasi-lD systems is crucial since it allows the identification of the 
frontier between the ballistic and the diffusive propagation of wavepackets in weakly disordered systems. Assuming 
that the transport regime remains coherent, a new class of scattering paths will yield an important contribution to 
the resis tance, known as the weak localization correction, which eventually turns the metallic state to an insulating 
one [ogl . llOd . llOll . Il02l |. The localization length ^ is the other physical length scale that defines such a transition, 
where the conductance decays exponentially with the system length as G '--^ Go exp(— L/^). 

Weak locali zation phenomena have been observed in multiwalled carbon nanotubes with diameter ranging from 
3 - 20 nm (Tol, ^M, ^QB, ^M,^M- Similarly, weak localization has been recently clearly unveiled in GNRs 
wit h widths i n the order of ~ 200 — 500 nm |l20l | . Additionally, transition to weak anti-localization has been reported 
[3, Il20l . Il2ll | as well. Weak anti-locaHzation (WAL) in graphene-based materials is argued to originate from some 
pseudospin-induced sign change of the quantum correction, similar to what is observed in systems with strong spin- 
orbit coupling (q^. [qsI. IM. [qbI . For the same kind of disorder potential, preservation of pseudospin symmetries might 
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jeopardize for GNRs with small width. We can expect that the effect of edge disorder and intrinsic defects (topological, 
vacancies, adsorbed impurities, ...) will play an increasing role as GNR width decays from ~ 20 nm down to '--^ 5 nm. 

It is thus genuinely important to evaluate the varying effects of disorder on quantum transport as the dimensionality 
or symmetries are changed. It is worth stressing that, in order to unveil weak localization effects, an external magnetic 
field is generally applied to tune the intensity of quantum corrections. Indeed, the phase of the quantum wavefunction 
is modified by the magnetic field (through a term giving the circulation of the potential vector along the scattering 
path), which reduces the probability of return to the origin and enhances the conductance (this phenomenon is 
known as the negative magnetoresistance effect). However, for low dimensional systems such as carbon nanotubes, 
it was demonstrated that magnetic field has also severe consequences on the ba.nd structure, so th at the resulting 
magnetofingerprints in locaHzation regimes become more complicated |l03l . 1 104 IiOTI . llOSl . Il09l . Ill0l | to follow. This 
will also apply to magnetoresistance effects in GNRs with width < 10 nm. 

Weak-disorder limit (W = 0.5) 
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FIG. 12: (a) Conductance for a single disorder configuration of a zigzag (solid blue line) and an armchair (dashed red line) 
GNR with width ~ 20 nm and for W = 0.5. Black lines correspond to ideal zigzag (solid line) and armchair (dashed lines) 
ribbons, (b) Configuration averaged (over ~ 400 samples) normalized conductance as a function of GNR length for both zigzag 
and armchair GNRs. The solid blue (dashed red) arrow shows the energy at which the calculations for the zGNR (aGNR) have 
been performed, (c) and (d) Same informations as for (a) and (b) but for a larger disorder strength {W = 2). Adapted from 
[Lherbier, A.; Biel, B.; Niquet, Y. M.; Roche, S. Phys. Rev. Lett. 2008, 100, 036803]. 

Let us illustrate the quantum localization effects in GNRs. As seen in Section II-B, zigzag type GNRs display very 
peculiar electronic properties, with wavefunctions sharply locaHzed along the r ibbon edges for energies in the vicinity 
of the CNP. Using a conventional Landauer-Biittiker approach (35l.[3^. [37l . ll22l |. one can explore the scaling properties 
of the quantum conductance of these systems. In Fig. fT2lfa) and Fig. fT2l^ b). the energy-dependent conductance for 
both zGNR and aGNR of width ^ 20 nm are shown for pure, weak disorder {W = 0.5) and strong disorder {W = 2) 
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limits. 

In the weak disorder limit {W = 0.5, Fig. [T2lfb)). aGNRs appear to be much less sensitive to disorder effects than 
zGNRs with the same width. As can be seen in Fig. fl^c) and Fig. [T2lfe). the averaged normalized conductances 
are exponentially damped, following (InG/Go) ~ L/^, where an average over ^ 400 different disorder configurations 
has been performed. From these calculations, ^ is found to be up to two orders of magnitude smaller in zigzag than 
in armchair ribbons in the low disorder limit {W = 0.5), for an energy value far from the close vicinity of the CNP 
(following jH, [IB, [s^l)- In contrast, for disorder strength as large as W = 2 (Fig. [T2lfd)). the localization lengths 
for both types of ribbons are almost equal, showing that edge symmetry does not play any role. This result can 
be understood by the lower transport dimensionality in the case of zigzag edge symmetry, driven by more confined 
wavefunctions [iB, [13, S H [SI, [sil . 

An important result o f mes oscopic physics is that there exists a fun damental relationship between £ei and ^ known 
as the Thouless relation [l23l |. In a strictly ID system, Thouless [l23l | demonstrated that 

e = 24i, (50) 
whereas for quasi-lD systems (with N±{E) conducting channels), the relation was generalized as 

^E) = mN±iE) - l)/2 + l]lei{E), (51) 

with (3 a factor dependent on the time- reversal symmetry [l24l |- Avriller and co-workers |l25l | recently confirmed 
numerically such relation by studying chemically doped metallic carbon nanotubes. 
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FIG. 13: Localization length of a (5,5) CNT. Dashed line: value obtained from Eqn. I|5ip by naively using the perturbative 
values for lei and N± based on the self-energy and DoS of the disorder-free system. Solid line: same value based on the 
non-perturbative quantities including disorder in the self-energy and DoS. The numerical values indicated by error bars are 
obtained from logarithmic sample averaging. Adapted from [Nemec, N., Richter, K.; Cuniberti, G. New J. Phys. 2008, 10, 
065014]. 

As it was shown in Ref. [ll4l |- Eqn. (|5T|) still holds to good precision near van Hove singularities on the condition 
that the non-perturbative effects of multiple scattering are correctly included in the calculation of the DoS and the 
number of channels iV_L as it was described before for the elastic mean free path. As can be seen in Fig. [I3l this 
correction greatly improves the agreement with the true value obtained via sample-averaging. The remaining deviation 
is predominantly caused by the fact that the different conduction channels have very different velocities and thereby 
very different elastic mean free paths, whereas the Eqn. fSTj) is based upon the assumption of equivalent channels 
with one common elastic mean free path. 

The expressions for the DoS and the localization length obtained and demonstrated before for carbon nanotubes are 
general enough to apply to graphene nanoribbons as well. The data in Fig. fTjF a) demonstrates that the values obtained 
analytically do indeed match the results of numerical sample averaging. For semiconducting CNTs and GNRs, the 
behavior near the gap is correctly reproduced, showing a slight disorder-induced reduction of the gap-width. 

As can be seen in Fig.fCTb). the expressions still hold for zigzag GNRs over most of the energy spectrum, and even 
qualitatively describe the effects of the edge state around the Fermi energy: the extremely high DoS of the Fermi 
energy is smeared out in the energy range and the localization length is drastically reduced not only in the narrow 
region of the low-dispersive edge state itself, but - due to the disorder-induced broadening in energy - in an extensive 
region around the Fermi energy. However, the expressions fail to describe the actual shape of the energy-broadened 
peak in the DoS, as well as the fianks of the suppressed region in the localization length. These strong deviations can 
be traced to the assumption of equivalent conduction channels that is made in Eqn. (|5T|) . Near the edge state, this 
assumption fails completely. 



D. Edge disorder in GNRs 



In contrast to two-dimensional graphene and carbon nanotubes, graphene nanoribbons are subject to chemical 
passivation and roughness at the edges. In some cases, the nature of the chemical groups that passivate the edges 
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FIG. 14: (a) Density of states g{E), number of channels A'^ch, elastic mean free path iei and localization length lioc of an 
armchair GNR of width Na = 20 under the influence of Anderson disorder {a^ = 0.2 eV). Dashed lines: g (E) and A'ch of the 
clean system. Solid lines: g {E) from Eq. II44II . A^ch obtained as transmission through cross section through infinite system with 
self-energy Eq. l|43|l and lengths obtained from these. Data with error bars: values obtained numerically by averaging over 
~ 180 samples of length 2000 to 20000^uc. (b) The same for a zigzag GNR of width Nz = 18. The most prominent feature is 
the disorder-broadened peak in the DoS caused by the edge state that causes a strong reduction in the elastic mean free path 
and the localization length. The analytical values fail to describe the true shape of this feature because the expressions neglect 
the difference in the various channels. Adapted from [Nemec, N., Richter, K.; Cuniberti, G. New J. Phys. 2008, 10, 065014]. 



(usually H) can be determined experimentally, thus enabling a certain control of the ribbon. Unfortunately, the 
spatial regularity of the edges is much more difficult to achieve and the state-of-the-art etching techniques cannot 
avoid roughness. Direct or indirect evidences of edge disorder have been observed in many and different graphene 
samples, independently of the technique exploited to fabricate them. 

In the literature, the edges of graphene sheets exfoliated by cleavage of highly oriented pyrolitic graphite have been 
directl y investigated by sc anning tunneling microscopy (STM) in association with s cann i ng tu nneling spectroscopy 
(STS) [HH [Ha, [Hi, [Sil and by direct contact atomic force microscopy (AFM) in [l30l. [l3ll| . Micro- Raman spec- 
troscop y ha s also proved to be a valuable tool for studying the armchair or zigzag orientation of the ribbon edges 
locally [l33 |. 

In all these experiments, the structure of the edges turned out to be very irregular, with alternation and mixing of 
zigzag and armchair terminations, protrusions and dents or more complex structures. In general, armchair segments 
are considerably longer than zigzag fragments, thus evidencing a lower stability for the latter. The measurement of the 
transport properties of lithographically etched ribbons (itI . [18| also provides indirect evidences of the edge roughness. 
In particular, hints of inactive edge regions and dependence of the maximum resistivity on the ribbons width have 
been related to a possible disorder on the edges. 

Similar conclusions are drawn for ultrathin epitaxial graphene grown on silicon carbide crystals (gI. Il33l |. In this 
case, an indication of edge disorder (due to roughness or chemisorbed molecules) comes from the lower-than-expected 
number of conductive modes, due to inactive edge regions, and from the behavior of weak localization. 

Innovative techniques for fabrication and etching of graphene ribbons have proved to reduce the edge irregularities 
considerably, besides allowing the realization of ultranarrow (few nm) structures. The chemical technique developed 
by Dai and co-workers [l^l enables the formation and selection of long and ultrasmooth nanoribbons. Despite the high 
spatial regularity, field-effect transistors based on these systems [l34l | have shown a non-negligible scattering r elated 
to disorder at edges. The very recent STM lithographic technique developed by Tapaszto and co-workers [l35l |. 
allows the patterning of ultranarrow structures with the possibility of choosing orientation and width with almost 
atomic precision. Nevertheless, STM measurements on few nm wide ribbons have revealed irregular oscillations in the 
electronic density of states, thus suggesting the possible presence of edge disorder and its importance for very narrow 
ribbons. 

The literature proposes several and different edge disorder models for graphene nanoribbons. The basic idea to start 
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FIG. 15: (a) Differential conductance of a metallic armchair ribbon (N=53) with a ~10% of vacancies on the edges. Couple 
of atoms have been removed on the most external rows over a length of about 1 fim. For the calculation, we adopted a simple 
tight-binding Hamiltonian with zero site energy and hopping energy 70 — —2.7 eV. (b) Differential conductance for the same 
system with Anderson disorder on the edges. The energy of the edge sites varies randomly in the range ±1/ with V =0.25, 0.5 
and 0.75 eV over a length of about 1 ^tm. 



from a tight-bind ing Hamiltonian for a clean arid reg ular ribbon (in general with hydrogen passivation on the edges), 
and then adding [30| or removing [ssl . Il36l . Il37l . Il38l | carbon atoms at t he edges, or varying the width of the system 
[l39l | to account for roughness, or introducing Anderson disorder |l40l |. The most proper way of adding/removing 
atoms at the edges avoids final configurations that might cause steric problems. In practice, only H-C-C-H groups can 
be removed or added at the edges of an armchair ribbon, while particular care must be adopted when the disorder 
goes deeper than the first row in zigzag ribbons. [Not all the authors take care of these constrains.] 

The transport pr oper ties of metallic aGNRs have been investigated in the presence of both roughness [35l . Il38l | and 
Anderson disorder |l4d | on the edges. In the case of vacancies on the edges, the differential conductance of the system 
is considerably reduced and even a weak disorder on the two external rows induces a localization process, see Fig. 
flSFa). The locaHzation length turns out to be particularly low around the neutrality point even for a 5% of vacancies. 
This is attributed to a small gap opening in correspondence of t he cr ossing of the two bands around E=0 [s^ and 
whose width is inversely proportional to the width of the ribbon [l38l |. When increasing the level of disorder on the 
two most external rows, the system can be envisioned as a sequence of metallic and semiconducting fragments and 
the localization length decreases considerably within the energy range that corresponds to the semiconductor gap, see 
Fig. [TBI Even outside this region, the conduction ability is seriously jeopardize d, w ith a localization length of few 
tens of nm. [Similar results have been also recently obtained by Evaldsson et al. |l37| .] 

The effect o f seq uences of large conducting and semiconductor armchair fragments has been investigated in detail by 
Martin et al. [139|. The length of each metallic fragment is such to preserve the band structure of the corresponding 
infinitely extended ribbon. From this perspective, it is possible to obtain an effective tight-binding Hamiltonian, 
where the eigenfunctions of the metaUic fragment play the role of "orbitals", with equispaced site energies, and the 
hopping energies between two subsequent fragments are related to the properties of the semiconductor fragment in 
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FIG. 16: One half the amplitude of the localization length for a metallic aGNR with A*" = 53 for a 5% (a) and a 50% (b) 
concentrations of defects at the edges. This figure is taken from [Areshkin, D. A.; Gunlycke, D.; White, C. T. Nana Lett. 2007, 
7, 204], by courtesy of C.T. White. 

between. The result is a one-dimensional impurity band insulator with conductivity 

cr sa g-V^-S^s/^ for T < T* (52) 

where Lav > W is the average length of the metallic grains. Eg is the energy gap of the semiconductor fragments, T 
is the temperature, T* « \jo\D/L1^ (with D the width of the ribbon) and a is a numerical coefficient of the order of 
1. 

In the case of weak Anderson disorder on the edges, the conductance of metallic aGNRs is only slightly affected, 
in particular within the energy region that corresponds to a single channel, see Fig. [rsTb). Thi s result can explained 
by considering the high kinetic energy of the states around the charge neutrality point [l4(]| |. which are therefore 
scarcely affected by potential fluctuations. Just to give an idea of the low sensitivity of these states to disorder, let us 
consider a constant potential on both edges. It turns out that the structure of the highest valence band and the lower 
conduction bands does not change around the charge neutrality point but for a small shift in energy and the rising 
of an energy gap whose width is always small and considerably suppressed for larger ribbons, see Fig. [iTl The linear 
energy dispersion is thus preserved and the gap is almost negligible because the states that correspond to the bands 
around the neutrality point are spread all over the section of the system and thus the effect of the edge potential is 
weaken by the averaging over the chains. 

The two types of disorder, always confined to the two outer rows of the metallic aGNRs, have different consequences. 
However, there is no inconsistency in this, since the nature of the perturbation is completely different in the two cases. 
Roughness tends to introduce semiconductor islands and then a gap. Weak Anderson disorder does not perturb the 
structure of the energy bands close to the neutrality point, therefore the backscattering is limited, at least for wide 
ribbons. 

In the case of roughness for semiconductor aGNRs, the system goes toward a strong localization re gime, as in 
the case of metallic ribbons. However, the disorder induces states in the region of the energy gap (s^. Il36l |. thus 
allowing current tunneling for not very long systems. Therefore, a compromise between the width and the length of 
the ribbon must be found in order to keep the semiconductor behavior of the strip. In particular, the ribbon must 
be short enough to preserve the extension of the states just outside the gap all over its length and long enough to 
prevent tunneling through the states induced within the energy gap. As re ported at the end of the section, this can 
considerably affect the lon/IoS ratio for GNRs-based field effect transistor |l4ll . Il42l . Il43l |. 

In the case of zigzag ribbons, the almost dispersionless edge states around the neutrality point play a key role 
in determining the properties of the disordered system. In the literature, different behaviors have been observed 
depending on the specific model of edge disorder, with results not always consistent. Areshkin and co-workers js^ 
have found that zGNRs are much less sensitive to edge disorder than aGNRs, especially in the energy region where 
only one conductive mode is active. For reasonably wide strip and for an erosion of 50% C atoms on the eight outer 
rows, the localization length turns out to be of the order of 10 /im (at least ten times the value for an aGNR with 
a much weaker disorder). To explain this, they consider that the states within the first conductive channel can be 
divided into two groups: almost dispersionless edge states within a very small energy range around the neutrality 
point E ^ and bulk states with energy outside this range. The width of the range is determined by the transverse 
width of the ribbons. The effect of roughness is mostly on the edge states, and this explains why the bulk states are 
not much affected by disorder. On the other hand, the conductance is expected to be much depressed around the 
neutrality point, where the effect of disorder is larger. This is only partially observed in the simulations js^, where the 
conductance seems anyway to keep a strong resistance to roughness with respect to the armchair case. This behavior 
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FIG. 17: Gap between the highest valence band and the lowest conduction band in the presence of a constant potential W on 
the two edges of the armchair ribbon for A'' =14,35,56,77,98,119,140. The larger is the strip, the smaller is the gap. 
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FIG. 18: Differential conductance of a zigzag ribbon (N=32) with Anderson disorder on the edges. The energy of the edge sites 
varies randomly in the range ±V with V =0.25, 0.5 and 0.75 eV over a length of about 1 ^m. 

is not understood and would deserve further investigation. 

A rather different conclusion is obtained by Querlioz and co-workers |l36l |. who evaluated and analyzed the wave- 
function and the density-of-states of zigzag ribbons. From these data the mobility edge is extracted. Their conclusion 
is that no ribbon has zero mobility edge and that the resistance to edge disorder for zGNRs vanishes as soon as more 
than one edge row is eroded. In this case, the roughness induces an Anderson insulator behavior independently of 
the zigzag or armchair orientation of the ribbon. This result, in agreement with other recent calculations [l38l |. seems 
to be in striking contrast with [sHi- More detailed analysis focusing on specific disorder models could unravel this 
ambiguity. 

Weak Anderson disorder on the edges leads to the opening of a gap in the conductance of zigzag ribbons [l4(]| |. 
Some results are reported in Fig. [181 A weak potential on the edges has a deep impact on the electronic spectrum and 
then on the transport properties around the CNF. The states corresponding to the two almost fiat bands at E ~ 
at the borders of the first Brillouin zone {k = zLn/a) are completely localized on the edges of the ribbon. Therefore, 
a potential on the edges moves (upward or downward) the energy of these states. As a consequence, the electrons 
are subject to backscattering, especially in the region where most of the onsite disorder energies are concentrated, i.e. 
around the neutrality point. This induces a gap in the conductance, whose width increases with the strength of the 
potential. 

We conclude this section by giving a look to the poss ible effects of edge roug hness on the behavior of graphene-based 
(field-effect or Schottky-barrier field effect) transistors [l4lLll42l . ll43l . fl43 . ll45l |. These systems are usually investigated 
with the help of Schrodinger-Foisson solvers, which allow the self-consistent evaluation of the electrostatic potential. 
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The results presented in the hterature agree that edge disorder reduces the Jon current and increases the loS current, 
thus making their ratio worse. The decreasing of Ion is due to quantum transport effects, as seen before, and to 
self-consistent electrostatic effects that compensate the vacancies with an accumulation of charge on the atoms close 
to it. The increasing of the leakage current is due to the tunneling of electrons through the states induced in the 
band gap region or through conducting zigzag fragments. The degrading of the /on/^off ratio is thus a serious issue 
for the efficiency of transistors based on very small ribbons. On th e other hand, larger ribbons wou ld re duce the gap 
tremendously. A recent theoretical paper by Ouyang et al. |l45l |. in agreement with experiments |l34| . shows that 
the presence of optical phonons (OPs) reduces the effect of edge irregularities considerably. In fact, electrons loose 
much energy by emitting an OP and they can hardly come back to the source after being reflected by the edges. As 
a consequence, elastic backscattering is only relevant when occurring close to the source, where electrons keep their 
energy before phonon emission. 



E. Minimum conductivity and charge mobility in graphene based systems 

In semiconducting materials (such as silicon nanowires [ll8l |). charge mobility is a very important quantity to assess 
the transport efflciency of the system and the corresponding device performances. By deflnition, the charge mobility 
reads: 

fi{E)=asc{E)/en{E), (53) 

where 

= e^p{E)v{E)£ei (54) 

is the semiclassical conductivity deduced from the Einstein formula, with p{E) the DoS, n{E) the charge density 
at energy E, i^i the elastic mean free path, and e the elementary charge. Close to the charge neutrality (or Dirac) 
point, the measured experimental conductivity of various samples was found to range within ~ 2 — 5e^/h although 
the charge mobility was changing by almost one order of magnitude [sl. IgI. ItI. I^. [qI. lid, fill, fl^. fisl. fl3|. This effect has 
been attributed to the change of charge density due to the doping from the substrate and/or contacts. 

On the theoretical side, the calculation of the Kubo conductivity for 2D graphene with short range disorder, and 
within the self-consistent Born approximation, yields a^^^ = 4e^/(7rft,) {h the Planck constant) for the two Dirac 
nodes [l46l |. which is typically I /tt srnalle r than all the experimental data. Numerical calculations using the Kubo 
formula confirm such prediction [l47l . Il48l | . Amazingly, as discussed in section HIE H this value also comes out in a 
completely ballistic transport regime, as a contact effect. 

Differently, by assuming that the Dirac fermion scat terin g is dominated by Coulomb scattering from ionized impuri- 
ties near the graphene plane, Nomura and MacDonald [l47l | could numerically reproduce the low energy dependence of 
the measured cha rge c onductivity, and found that cr™" ~ e'^/h close to the Dirac point, in better agreement with most 
experiments. In |l47l |. the authors used the full quantum approach of the Kubo formula, describing the long-range 
disorder effects via a screened Coulomb potential, and performing a finite-size scaling analysis. 

Other calculations within the semi-classical Boltzmann approach or th e Landauer approach have also reported on 
the effect of screened Coulomb potential on charge conductivity [13, Il49l | with similar conclusions, and although the 
existence of a true universal minimum conductivity could not be rigorously answered, an interpretation in terms 
of saturat ion o f the conductivity due to charged impurity induced inhomogeneities, occurring at low densities, was 
proposed |l49l |. 

These results suggest that the intrinsic disorder in graphene could be of electrostatic (Coulomb scattering) nature, 
likely due to charges trapped in the oxide. One also notes that the effect of graphene plane deformation modes (known 
as ripples), frozen when the exfoliated la yer i s dep osited onto an oxide layer, has been investigated by introducing 
an effective random gauge potential jo^, Il50l . Il5l| . The authors found that a temperature independent minimal 
conductivity will take place, with full suppression of localization effects owing to the absence of intervalley scattering 
processes, in agreement with temperature-dependent experiments (20| . 

Concerning the contribution of quantum interference effects and the transition to an Anderson type localized regime, 
the issue is still controversial. In presence of strong intervalley scattering processes, which is best reaHzed for short 
range disorder potential, conventional weak locaHzation phenomena have been predicted [92. ,152i | and observed [12o| . 
The preservation of the pse udospin s ymmetry also allows for the manifestation of a spectacular symmetry-dependent 
anti- localization effect (9^. Il2(]| . Il52l |. The transition to a localized Anderson type regime is more debated |9^. [97|. 
with to date no experimental evidences in such material. The question is how quantum interferences effects and 
localization phenomena for Dirac fermions depend on t he u nderlying disorder potential characteristics, and whether 
the conventional 2D scaling theory of localization |lOll . llod | is applicable or not in today's graphene materials. 
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Several authors have challenged the single parameter scaling theory of Anderson localization by studying the so- 
called beta function f]{g) = d\ng/d\nL, with g the dimensionless conductance and L the system size. By computing 
the scaling behavior of /3{g), the locaHzation versus delocalization nature of electronic states can be analyzed. In 
conventional two-dimensional systems, the theory predicts that al l sta t es a. re localized, independently of the disorder 
characteristics and provided time reversal symmetry is preserved [lOlL llod |. 

Recent numerical studies claim that in presence of short range diso rder (Anderson-type potential) and intervalley 
scattering, all states are indeed locaHzed even for Dirac fermions [l48l |. However, it is interesting to note the typical 
values obtained for transport length scales. By varying the disorder strength from W ~ 3.4eV to ~ 2eV, the 
elastic mean free path was found to range from a few nm to several tens of nm close to the Dirac point fFigfTOl-left 
panel (b)), wh ereas the locaHzation length £, given by ^ = £eexp{TTasc/Go) upscaled from 20nm to 10//m (not shown 
here, see Ref. [l48l |). The energy-dependence of ^ was also shown to be strongly driven by that of the semi classical 
conductivity with a minimum localization length at the Dirac point |l48l |. 

In case of long range disorder, the situation is more complex. Indeed, in absence of intervalley scattering, Dirac 
fermions cannot be trapped by a potential well, irrespective of the well depth. This suggests the robustness of states 
against an i nsula ting tendency. A different kind of scaHng behavior of the conductance at the Dirac point was first 
proposed in |l53l |. The key result was the occurrence of a quantum critical point at half filling giving some universal 
value of the conductivity of the order of e'^/h, so in contradiction with a localized nature giving a zero conductivity in 
the thermodynamic limit. Other numerical studies (98l . ll54| found a different scaling flow for the beta function of the 
Dirac model, thus indicating that all states remain delocaHzed whatever the strength of the underlying disorder is, 
but with a conductivity upscaling with length up to infinity. The situati on c hanges in the presence of ripples, which 
can be described by random vector potentials. Nomura and co-workers |l5Cll | further deepen the scaling behavior of 
both transport coefficients <Jxx and cr^y and found that massless Dirac fermions will exhibit a critical behavior similar 
to that of the quantum Hall transition point, but in absence of uniform magnetic field. 

An important observation is that to date, the theoretical description of disorder in graphene layers (either deposited 
on a substrate or suspended in between contacts) is mostly achieved at a phenomenological level. Usually, short 
range and long range scattering potential are described by some onsite potential fluctuations (Anderson-type) and 
by Gaussian correlated potentials, respectively. Alt hough the study of deformation modes known as ripples is at 
the origin of peculiar transport predictions [Oa . Il55l |. a realistic description of Coulomb scatters is howev er ne eded 
to allow a true experimental exploration of locaHzation phenomena in massless Dirac fermions. Novikov |l56l | flrst 
discussed about the possible asymmetry in the transport cross section for a Dirac electron scattering off a positively or 
negatively charged Coulo mb i mpurity. Self consistent RPA-Boltzmann theory also found Coulomb scattering induced 
conductivity asymmetry [l57j . without the possibility however to tackle with localization effects. A recent theoretical 
study has investigated quantum coherent transport and transport scaling lengths for intentionaHy chemically doped 
(and disordered) 2D graphene layers with a reaHstic and self consistent description of impurity scattering potentials 
|158l |. By incorporating substitutional boron (or nitrogen) impurities, elastic mean free paths, as well as the semi- 
classical conductivity and charge mobilities were numerically estimated by the Kubo approach. Some onset of quantum 
interference effects was also observed, even at the Dirac point, but this contribution was found too small to explore 
the possible underlying scaling behavior of the beta function, even in situation of strong doping (such as 4%). 

Finally, although conductivity is a well defined quantity, 2D graphene manifests specific properties that make the 
use of formula Eq. (|53| somehow ill-defined. Indeed, when the energy of charge carriers approaches the Dirac point, 
the semiclassical conductivity remains finite, whereas /i ^ oo, as the charge density n — > 0. 

In Fig. [191 show the result of a numerical calculation using the Kubo approach (see [l48l | for details). The 
evolution of ^ is shown as a function of the Anderson disorder strength W. The energy-dependence of ^J,{E) and 
£ei{E) are found to be similar. In the close vicinity of the Dirac point, the downscaling of fi with W fol lows the 
Fermi Golden Rule prediction, while it diverges when approaching the Dirac point. Experimental data from |l59l | are 
also shown in Fig. [19] (right) for comparison. Different samples with varying quality show a similar trend, although 
patterns from different samples can substantially differ in shape indicating fiuctuations in the disorder characteristics. 

It is clear that by further reducing the dimensionality of the graphene material, charge mobilities wiH be reduced. 
Recent transpo rt ni easurements on nanoribbon-based field effect transistor show mobilities in order of ^ ^ 100 — 
SOOcm^V^s-i [Isl, which are thus reduced in comparison with the 2D graphene measurements (reported values can 
be as large as a few lOCOOOcm^V^s-^Q, S, i, [13, [13, [H, Q)- However, the lateral size reduction aHows for a 
larger energy gap, which ensures more efficient field effect efficiency. 

IV. CONCLUSION 

To conclude, in this review we have reported on the basics of electronic and transport properties in low dimensional 
carbon-based materials including 2D graphene, graphene nanoribbons and carbon nanotubes. It has been shown that 
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FIG. 19: Left: charge mobility as a function of carrier density (a) and mean free path (b) for W = 1, 1.5,2,2.5 (from top to 
bottom), from [Lherbier, A.; Persson, M. P.; Niquet, Y. M.; Triozon, F.; Roche, S. Phys. Rev. B. 2008, 77, 085301]. Right: 
experimental mobility from [Tan, Y. W.; Zhang, Y.; Bolotin, K.; Zhao, Y.; Adam, S.; Hwang, E. H.; Das Sarma, S.; Stormer, 
H. L.; Kim, P. Phys. Rev. Lett. 2007, 99, 246803], by courtesy of Philip Kim. 

although nanotubes and nanoribbons share similar electronic confinement properties due to their nanoscale lateral 
sizes, the effects of boundary conditions in the perpendicular direction with respect to the system axis trigger very 
different transport features when disorder is included. Close to the charge neutrality point, the robustness of armchair 
metallic nanotubes against disorder is absent in nanoribbons, which cannot be classified in the family of ID massless 
fermions, owing to edge- induced gap openings. Nanoribbons with zig-zag symmetries are even more spectacularly 
sensitive to disorder owing to the edge-states driven lower transport dimensionality. In contrast, for charge carrier 
energies lying in the higher energy subbands, the properties of nanotubes and ribbons present similar features, with 
strong energy dependence of elastic mean free paths and localization phenomena. 

Additionally, the transition from a quasi-lD to a true 2D system results in strong damping of disorder effects, with 
enhanced elastic mean free paths together with strong damping of quantum interferences. In particular, the study 
with Anderson disorder demonstrates that even in the strongest case of short range scattering potential (with possible 
short range potential fiuctuations as large as leV) , the computed 2D localization lengths remain in the range of several 
hundreds nanometers to microns. One can thus conclude that to observe weak and strong localization regimes, the 
presence of edges as well as a reduced lateral size are essential factors. 

Finally, the possibility to produce and control defect densities either through intentional doping or by irradiation 
techniques (that produce vacancies-type defe cts') could open spectacular avenues to explore quantum transport phe- 
nomena (including quantum Hall effects [l6(]l |) in low dimensional materials, for which a realistic description of both 
underlying electronic structures as well as superimposed disorder potentials would be theoretically possible, allowing 
unprecedented exploration of experimental data at a quantitative level. 
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